Relativistic Brownian motion: Prom a microscopic binary collision model to the 

Langevin equation 
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The Langevin equation (LE) for the one-dimensional relativistic Brownian motion is derived 
from a microscopic collision model. The model assumes that a heavy point-like Brownian particle 
interacts with the lighter heat bath particles via elastic hard-core collisions. First, the commonly 
known, non-relativistic LE is deduced from this model, by taking into account the non-relativistic 
conservation laws for momentum and kinetic energy. Subsequently, this procedure is generalized 
to the relativistic case. There, it is found that the relativistic stochastic force is still (5-correlated 
(white noise) but does no longer correspond to a Gaussian white noise process. Explicit results for 
the friction and momentum-space diffusion coefficients are presented and discussed. 

PACS numbers: 02.50.Ey, 05.40.-a, 05.40.Jc, 47.75.-|-f 



I. INTRODUCTION 



The theories of the non-relativistic Brownian motion 
and special relativity were introduced more than 100 
years ago H S H S S IS H- Since then, they have 
become cornerstones for our understan ding of a wide 
range of physical processes 0, 0, 0, E|. This fact 
notwithstanding, the unification of both concepts poses 
a theoretical chall eng e still nowadays (classical references 
are [13, |3 11 . fla. [l ^ ITsIl: recent contributions in- 
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potential applications in high-energy physics and astro- 
physics are considered in |33, 0, Issl l3a. l3^ Issf). The 

relatively slow progress in this field can be attributed 
to the severe difficulties that arise when one tries to 
describe iV-body systems in a relativistically consistent 
manner (3£, ^40, 41, 43j 03- Due to this reason, the 
derivation of relativistic Langevin equations (LEs) from 
an underlying microscopic model has remained an un- 
solved issue until now j67| . However, in the present paper 
we aim to provide a solution to this problem. 

More precisely, by considering quasi-elastic, binary 
collisions between the Brownian and heat bath parti- 
cles 0, we are able to treat the heat bath in 
a fully relativistic manner without having to account 
for the exact details of the relativistic iV-body interac- 
tions. As shown in Sec. ^ for a non-relativistic frame- 
work this approach yields the well-known non-relativistic 
LE with Gaussian white noise as well as the correct 
fluctuation-dissipation theorem (the Einstein-Sutherland 
relation Q]). In Sec. lIIII the method is transferred to the 
relativistic case, leading to the main result of this pa- 
per, the relativistic LE H32I) . Remarkably, the relativistic 
stochastic force is also (^-correlated ('white') but no longer 
of Gaussian (or Wiener ^^) type. Compared with the 
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non-relativistic Brownian motion, this is the most impor- 
tant difference. Furthermore, we obtain explicit represen- 
tations of friction and (momentum)-diffusion coefficients 
in terms of expectation values with respect to the heat 
bath distribution (see also Add. 



II. NON-RELATIVISTIC BROWNIAN 
MOTIONS 



The objective of this section is to recover the well- 
known non-relativistic LEs from a simple microscopic col- 
lision model for Brownian motions. As shown by several 
authors in the past 4&, ^ , non-relativistic 

LEs can also be derived by considering a bath of har- 
monic oscillators with canonical phase space distribu- 
tion. Unfortunately, it is problematic to transfer this ap- 
proach to the relativistic case, because any instantaneous 
linear (or nonlinear) interactions between Brownian and 
heat particles would violate the basic principles of spe- 
cial relativity "JT, '42] . To circumvent this problem, we 
will pursue a different method here, using only the (non- 
)relativistic microscopic conservation laws for energy and 
momentum, respectively, known to hold for elastic point- 
like, binary collisions (contact interactions 42]). ConcCT 
tually, our approach is related to that of Pechukas 
and Pechukas-Tsonchev [45| . who considered a similar 
model in the context of non-relativistic quantum Brow- 
nian motion [s^ . Analogous approaches are also known 
from unimolecular rate theory, see e.g. Sec. V in [s^ . 



A. Microscopic model 

For the sake of simplicity only, we will restrict ourselves 
throughout to the one-dimensional (Id) case. Generaliza- 
tions to higher space dimensions are in principle straight- 
forward, but certain calculations will become much more 
cumbersome (cf. comments at the end of App. 0. To 



2 



start out, consider the following situation in the inertial 
laboratory frame Eq- A large one-dimensional box vol- 
ume V = [— i/2,L/2] contains an ideal non-relativistic 
gas, consisting of N small point-like particles with iden- 
tical masses m. The gas particles - referred to as 'heat 
bath' hereafter - surround a Brownian particle of mass 
M ^ m. Due to frequent elastic collisions with heat 
bath particles, the Brownian particle performs stochas- 
tic motions. 



1. Heat bath 

The coordinates and momenta of the heat bath parti- 
cles are denoted by Xr G [— L/2, L/2\ and pr £ (—00; 00), 
respectively, where r = 1,...,N. As usual, we make 
the following simplifying assumption concerning the heat 
bath: The probability density function (PDF) of the heat 
bath particles is a spatially homogeneous Maxwell distri- 
bution, i.e., at each time t > 0, the PDF reads 



N N 



Pr 



2mkT / ' 



(1) 



where k is the Boltzmann constant, T the temperature, 
and A = (27rmfcr)~^/^. Thus, it is implicitly assumed 
that: 

• the heat bath particles are independently and iden- 
tically distributed; 

• the collisions with the Brownian particle do not sig- 
nificantly alter the bath distribution. 

These assumptions are justified, if the collisions between 
the gas particles rapidly reestablish a spatially homoge- 
neous bath distribution. 



B. Derivation of the Langevin equation 

The total momentum change 5P of the Brownian par- 
ticle within the time interval r can be written as 

N 

6P{t) = P{t + r) - P{t) = ^r{t, t), (5) 



where Ir{t, r) G {0, 1} is the indicator function for a col- 
lision with the heat bath particle r during the interval 
[t,t + T]; i.e. Ir{t,T) = 1 if a coUision has occurred, and, 
otherwise, Ir{t, r) = 0. In the Id case, the collision indi- 
cator can be written explicitly as 

Ir{t,T) = e{X ~ Xr) Oix'^ - X') e{Vr - V) + 

e{xr - X) e{X' -x'^) OiV -Vr), (6a) 



where X = X (t) , Xr = Xr{t), and 
X' = X + Vt, 



(6b) 



are the projected particle positions at time t + t. The 
Heaviside-function is defined by 



e(x) 



The expectation of the collision indicator with respect 
to the bath distribution, denoted by (/^(i, r))j^, gives 
the probability that the bath particle r collides with the 
Brownian particle between t and t + t. As shown in 
App.0 in the limit r ^ 0, one finds 




(7a) 



2. Kinematics of single collision events 

The momentum and energy balance per (elastic) colli- 
sion reads 



E 



E 



P+p = P + p. 



(2) 



Here and below, capital letters refer to the Brownian par- 
ticle and small letters to particles forming the heat bath; 
quantities without (with) hat-symbols refer to the state 
before (after) the collision. In the non-relativistic case, 
we have, e.g., before the collision 



MV, 



p = mv, 



E 



2M' 



2m' 



(3) 



where v and V denote the velocities. Taking into account 
both conservation of momentum and (kinetic) energy, 
one finds that the change AP = P — P oi the Brown- 
ian particle's momentum per single collision is given by 



AP: 



-2m 



M + m 



P 



2M 

M + m 



(4) 



with function C{P) = C(V{P)) given by the integral 
formula 



C{V) 



1 



dVr {Vr - V) fl(Vr) + 
V 

AVr {V ~ Vr) jl{Vr). (7b) 



Here, fi,{vr) is the one-particle velocity PDF of a heat 
bath particle. We anticipate that Eqs. © and O remain 
valid in the relativistic case as well, but then one has to 
insert the relativistic bath distribution in Eq. 17bp . 

However, in order to recover from Eqs. ©-{7)1 the well- 
known non-relativistic LE, we still have to make a num- 
ber of simplifying assumptions: 

(i) The time interval r is sufficiently small, so that 
\SP/P\ <C 1. In particular, r is supposed to be so small 
that there occurs at most only one collision between the 
Brownian particle and a specific heat bath particle r. 
One the other hand, the time interval r should still be 
large enough, so that the total number of collisions within 
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T is larger than 1. These requirements can be fulfilled 
simultaneously only if m/M <^ 1. 

(ii) Collisions occurring within [t, t + r] can be viewed 
as independent events. 

(iii) Finally, we will (have to) assume that 

{[pr Irit.rW)^ - {jP, Ir{t,T))^ 

= (pi), CiP) J (8) 



for j = 1,2,.... Given the explicit representation of the 
indicator function (|6a|l . it is in principle straightforward 
to check the quality of the approximation (jHJ, if a bath 
distribution has been specified. 

As we shall see immediately, the assumptions (i)-(iii) 
are necessary and sufficient for deriving the well-known 
non-relativistic LE from Eqs. ©-0. Upon inserting 
Eq. I0J into |(SJ) and dividing by t we find 



5P{t) 



2m 



T ^ m + M 

r—l 



Irit,T) 



P 



1 ^ 



2M 



T ^ M + TO 

r— 1 



Pr Ir{t,T). 



(9) 



The first term on the rhs. in Eq. ^ can be identified as 
the 'friction' term, whereas the second term represents 
'noise'. On the rhs. of Eq. it was assumed that for 
each collision occurring within [t,t + t], the 'initial' mo- 
mentum of the Brownian particle is approximately equal 
to some suitably chosen value P{t') with t' e [t,t + T], cf. 
the assumption (i) above and the discussion at the end 
of this section. 

The next step en route to the conventional LE consists 
in replacing the square bracket expression in Eq. (|5J| by 
the averaged friction coefficient 



1 ^ 



2m 



M 



(10a) 



Since it was assumed that the heat bath particles are 
independently and identically distributed, we can rewrite 
this as 

, , N 2m , , , , 

for some r g {1, . . . , N}. The coefficient i/q can be in- 
terpreted as an average collision rate weighted by some 
mass ratio. Inserting Eq. (jTaf) into Eq. Ijl0b|l yields 



2to 



m + M 



C{P), 



(11a) 



where rib = N/L is the density of the bath particles. In 
the case of the Maxwell distribution, we can evaluate the 



integral (|7b|l . and find 

/ kT V'^ 

P 



m f P 



2M 



erf 



2kT \M 

1/2 P 



/ TO r 
\2kf) M 



(lib) 



In particular, setting (see App. ^ 



(12) 



corresponds to the commonly used Stokes approximation. 
It then remains to analyze the 'noise force' 



1 ^ 



2M 



T ^ — ' M + m 

r—l 



(13) 



corresponding to the last term in Eq. jnj. The momen- 
tum dependence of the noise enters through the implicit 
P-dependence of the collision indicator functions /^(i, t). 
To keep subsequent formulae as compact as possible, we 
shall use the abbreviation ^(t) = ^{P,t) in the remain- 
der. Then, averaging over the bath distribution and 
using Eqs. lO, we find for the mean value 

(e(i))b = 0. (14a) 

Furthermore, assuming mutual independence of the col- 
lisions, the correlation function is obtained as 

2 N 



\ M + TO 

m Sts^ f 2M 
~ T'^\M + m 

unu 5ts ( 2m2 



2 N 



J2 {Ir{t,T))^ 



T \M + TO 



(14b) 



with 6ts G {0, 1} denoting the Kronecker-symbol. To ob- 
tain the second line, we have used that l'^{t, t) = Ir{t, r), 
and the simplifying assumption ||SJ) that Ir{t,T) and Pr 
are (approximately) independent random variables with 
respect to the bath distribution. 

Similar to Eq. Ill4bl) . also the higher correlation func- 
tions are determined by the corresponding moments of 
the Gaussian marginal bath distribution Thus, un- 
der the above assumptions (i)-(iii), the non-relativistic 
stochastic force ^ {t) corresponds to Gaussian white noise 
(or a Wiener process respectively). 

Finally, by substituting vq from Eq (|11|1 for the square 
bracket expression in Eq. ^ and formally letting r — > 0, 
we recover from Eq. ^ the well-known non-relativistic 

LE lEllllii 



p = -MP)p+m, 



(15a) 



4 



where £,{t) = ^{P,t) is a momentum-dependent Gaussian 
white noise force, characterized by 



with (momentum-space) diffusion coefficient 

M2 



Do{P) 



M + m 



MP) kT. 



(15b) 
(15c) 



(15d) 
s) for 



To obtain Eq. H15c|) . we used that Sgt/r —^5{t 
r — > 0, where S{t — s) is the Dirac-function. 

In the fimit m/M 0, Eq. (|15d|) reduces to the 
standard fluctuation- dissipation theorem Dq — Mv^kT 
56]. However, vq and Dq are constants only if one 
adopts the Stokes approximation ((T^ . cf. App.^ If one 
goes beyond the Stokes approximation, then the noise 
in Eqs. (|15|l becomes muhiphcative with respect to P, 
and, therefore, Eqs. l(TK)l must be complemented by a 
discretization rule in this case Is^ Is^ 1581 Is^ Isol l6ll 
m m ill. As discussed in [slHealealel^ only for 
the post-point discretization rule, corresponding to the 
choice vo{P) = vo{P{t + r)) and Dq{P) = Do{P{t + r)) 
on the rhs. of Eq. H15al) . one recovers the Maxwellian 
PDF 



$oo(P) 



1 



2-KMkT 



1/2 



cxp 



2MkT 



(16) 



as the stationary momentum distribution of the Brow- 
nian particle in the limit t ^ oo (assuming that 
m/M 0). 



III. RELATIVISTIC BROWNIAN MOTIONS 

We shall now apply an analogous reasoning to obtain 
a relativistic LE. For this purpose we consider an iner- 
tial (laboratory) frame Eq with time coordinate t, as e.g. 
measured by an atomic clock resting in Eq. 



A. Microscopic model 

The basic constituents of the microscopic model are the 
same as those outlined in Sec. Ill Al but in addition we 
now have to consider a relativistic heat bath distribution 
and must consistently take into account the relativistic 
collision kinematics. 



1. Relativistic tieat bath, 

In the relativistic case, we postulate analogous to 
Eq. that, with respect to Sq, the heat bath distri- 
bution is stationary, spatially homogeneous, and inde- 
pendent, so that the PDF can be written in the product 



form 



N 



,Pn) = L'^Ylf^Pr 



(17a) 



r=l 



As marginal one-particle momentum PDFs, we will 
now consider the ?7-generalized Jiittner-Maxwell distri- 
butions 113, 113, reading: 



e{p)-n 



exp 



kT 



?7>0, 



(17b) 



where p G (— oo,-|-oo), and e{p) denotes the relativistic 
kinetic energy of a heat bath particle. The normalization 
constant Mri is determined by the condition 



1 = 



dp/b'b). 



(17c) 



For -q — 0, Eq. (|17b|l reduces to the standard Jiittner- 
Maxwell distribution "6^. On the other hand, as dis- 
cussed recently ^] , the PDF with rj = 1 appears to 
be conserved in relativistic elastic binary collisions. In 
general, however, the arguments and results presented 
below remain valid for arbitrary one-particle momentum 
PDFs /b(p), i.e., also for momentum distributions other 
than the 77-generalized Jiittner PDFs Ijl7b|l . 



2. Relativistic collision kinematics 

Using natural units such that c = 1, relativistic kinetic 
energy, momentum and velocity are related by 

1 /2 

p — mvj{v), e(p)=(TO^+p^) , (18a) 

P^MVj{V), E{P)^{M^ + P^Y''^ , (18b) 

— 1/2 

where 7(1;) = (l — v^) . As before, capital letters re- 
fer to the Brownian particle. Inserting Eqs. H18|) into the 
conservation laws Q, and solving for P, one finds [IJ 



2uE~ {l + u^)P 



where 



u{p,P) 



E 



(19) 



(20) 



is the center-of-mass velocity. Hence, the momentum 
change AP = P — P of the Brownian particle in a single 
collision is given by 



AP = -- 



1 - ii2 E + e 



P 



E 



1 



(21) 



In the non-relativistic limit case, where <C 1, — M 
and e ~ TO, this reduces to Eq. 
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B. Derivation of the Langevin equation 



Inserting Eq. (|^ into Eq. 
tivistic analogon of Eq. (O as 



6P{t) 



1 ^ 

-T- 



ISl, one obtains the rela- 



P 



r=l 
N 



-E- 

T ^ 1 



-Ul E + Er 
2 E 



T ^\-ui. E + Er 
r—1 ' 



Ir(t,T) 

PrIr{t,T), (22) 



where = u{pr,P) and = £{pr)- Formally, the colli- 
sion indicator Ir{t, r) is still determined by Eqs. © and 
(UJ, but differences arise due to the fact that we have to 
use V = P/{M'^ + p2)i/2 a,nd a relativistic bath distri- 
bution now. 

Analogous to the non-relativistic case, we can identify 
the first term on the rhs. of Eq. (|22|l as friction, and 
introduce an averaged friction coefhcient by 



1 ^ / 



N / 2 



ui E + e, 



Ir{t,T) 



l-ui E + e, 



Ir{t,T) 



(23) 



for some r G {1, . . . , N}. Next, applying a product ap- 
proximation similar to (O, we obtain 



y{P) 



N 

T 



□ 



C{P) 



2 er 



l-ui E 



(24) 



where nb = N/L is the density of the heat bath par- 
ticles, and C{P) is determined by Eq. I)7b|l . Figure 
shows the P-dependence of v{P)/[nbC{P)] for the bath 
distributions from Eq. (|17b|l . This momentum depen- 
dence is induced by the appearance of Ur = u{pr, P) 
and E — E{P) in the expectation value on the rhs. of 
Eq. H24I) . Furthermore, the shape of the one-particle col- 
lision coefficient C{P) — {Ir{t,T))^L/T is depicted in 
Fig. [3 As one would intuitively expect, the friction co- 
efficient grows with the temperature T of the heat bath 
(at constant P) as well as with the absolute momentum 
of the Brownian particle (at constant T). In the non- 
relativistic limit case, where <C 1, i? ~ Af and e 2± m, 
the relativistic friction coefficient i'(P) from Eq. 1)24(1 re- 
duces to the non-relativistic result i^oiP) from Eq. 

At this point, it might be worthwhile to emphasize 
once again that product approximations of the form 

{GiXr,Pr) Ir{t,T))^ ~ {G{Xr,Pr))^ {Ir{t,T))^, (25) 

as employed in Eq. © and also in the first line of 
Eq. (|24|l . can in principle be omitted by using the ex- 
plicit representation of the collision indicator and the 
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FIG. 1: The momentum-dependent, relativistic friction 
coefficient v{P), divided by the total mean collision rate, 
!/(P)/[nbC(P)], as calculated numerically for two different 
heat bath distributions f^(p) and two different bath temper- 
atures, is depicted versus the scaled momentum P. The solid 
lines refer to the standard Jiittner distribution with r] = 0, 
and the dotted lines to 77 = 1 in Eq. llTbl . (a) Weakly rela- 
tivistic heat hath. In the limit kT <^ mc? the bath distribu- 
tions l)17bfl approach a Maxwellian, and therefore the results 
for different r] practically coincide. In particular, for P = 
the non-relativistic result is recovered, (b) Strongly relativistic 
heat bath. The friction coefficient increases with the temper- 
ature of the heat bath. 



Eqs. HA6|I of the Appendix; if one opts to avoid such ap- 
proximations then the accuracy of the Langevin model 
increases (note that this statement applies to the non- 
relativistic case, too). However, in the following we shall 
continue to use Eq. H25|l in order to obtain a relativistic 
LE that is on an equal footing with the non-relativistic 

LE (nni. 

For this purpose, we interpret the second term on the 
rhs. of Eq. H22|) as 'noise', defining 



1 ^ 



E 



PrIr{t,T). (26) 



Averaging over the bath distribution f^ , one finds for 
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FIG. 2: Relativistic one-particle collision coefficient C{P) — 
(/r(t, t))j^ L/r, numerically calculated for the same parame- 
ters as in Fig. The solid lines refer to a standard Jiittner 
distribution with rj = 0, and the dotted lines to 77 = 1 in 
Eq. Ill7bll . (a) Weakly relativistic heat bath. At small tem- 
peratures, the zero- value C(0) is approximately equal to the 
non-relativistic Stokes value \/ kT/ (27rm). (b) Strongly rela- 
tivistic heat bath. For \P\ 00 the coefficient C{P) converges 
to 1/2. 



FIG. 3: Mean value of the relativistic stochastic force, 
fi{P) = (x(^))b' calculated numerically for two different heat 
bath distributions (p) and two different bath temperatures. 
Solid lines refer to a standard Jiittner distribution with 7; = 0, 
and dotted lines to 7; = 1 in Eq. Ill7bll . 



Let us next take a closer look at the covariance function 



the mean value 



M(P) ^ (x(P,t))b 

N / 2 E r , \ 

= — (-j J TT- PrIr[t,T)) 

T \l - uf. E + er I ^ 

T \l - uil: E + er I ^ 
^ n^C{P) /-^^—p \ . (27) 

In contrast to the non-relativistic case, the mean value 
/i of the relativistic Langevin force x(€) = x{Pj t) 
depends on the momentum P of the Brownian par- 
ticle. This can be attributed to the appearance of 
ul = {P + p^)"^ / {E + er)^ in Eq. As shown in 

Fig. 01 the quantity n{P) / [ubC (P)] is positive for P > 
and negative for P < 0. Thus, on average, the relativistic 
stochastic force tends to accelerate particles in the direc- 
tion of their motion, but this effect is compensated by 
the increase of the friction coefficient ^{P) at high values 

of F, cf. Fig. n 



a,s ^ {[X{t)~{x{t)k] [x(^)-(x(s)),])b- (28) 



In the non-relativistic case, the stochastic force possesses 
a vanishing mean value, (C(O)b ~ According to 
Eq. (|27|l . this is no longer the case for the relativistic 
noise force, (x(^))b 7^ 0- order to explicitly calculate 
ats for the relativistic case, it is convenient to introduce 
the abbreviation 



2 E 



Assuming, as before, that collisions can be viewed as in- 
dependent events, the correlation function H28|) vanishes 



7 



at non-equal times t ^ s, and we thus find 



E3 ^ 

Cts — Ots 



123 5 - 



- KrIr{t,T) 
T ^ — ' 

r=l 

^ r— 1 

N N 



St. 



N 



~ ^2 



rE («')b(^(^.-))b- 

(/,(t,T))bK)bE('^J-)b(^(^'^)) 

^ r t2 

^A^(^)E^M^)}- 



(29) 



From this, we obtain 



T T 

Ei Sts 



N 

5ts 2/ 



TV 



(30) 



The last term vanishes, if we consider the thermodynamic 
limit (TDL) of an infinite heat bath, i.e., TV, L — > cx) such 
that rib = N/ L =:constant. Thus, reinserting the explicit 
expression for k^, we obtain in this limit 



E 



r \\\ — ui L + Er 



(31) 



In principle, any higher correlation function can be cal- 
culated in the same manner. It is also evident that the 
noise force is non-Gaussian, because the relativistic bath 
distribution fh{Pr) that determines the averages ( • )b ~ 
and, thus, the noise correlations ~ is non-Gaussian. 
Finally, by substituting the averaged friction coefficient 



(32a) 



for the square bracket term in Eq. (|22|) . imposing the 
TDL for the bath and letting t in Eq. (jUl), we 
obtain the relativistic LE 



P=-v(P)P + x{t), 



(32b) 



where, in view of approximation H25|l . the non-Gaussian 
momentum-dependent noise force xi^) = x{Pt^) is char- 



acterized by the mean 

= "b C{P) 
and the covariance 



E 



l-ui E + er 



(32c) 



a{t,s) ^ ([x(t)-(x(i))b] [x(^)-(x(.'^))b])b 

= 2D(P) 5{t - s), (32d) 

with the (momentum-space) diffusion coefficient given by 

2\ 



Pi 



(32e) 



In Fig.Elthe ratio P'(P)/[nbC(P)] is plotted for the same 
parameters as in Figs. ^ and |3| As it is evident from the 
diagrams, this quantity increases with temperature T and 
absolute momentum P of the Brownian particle. 
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FIG. 4: Relativistic diffusion coefRcients D{P) calculated 
numerically for different heat bath distributions /b(p) for (a) 
weakly and (b) strongly relativistic heat. The solid lines refer 
to a standard Jiittner distribution with — 0, and the dotted 
lines to ?7 = 1 in Eq. (tTTbl . 



IV. RESUME 

We conclude the derivation of the relativistic LE with 
a set of general remarks: 

(i) While deriving the relativistic LE (|32|l . we made 
use of the stationarity, independence, and homogeneity 
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of the bath distribution (|17a|l : we did not, however, rely 
on the specific properties of the marginal momentum 
PDF. Hence, the above results hold true for arbitrary 
one-particle momentum distributions f^ip)- 

(ii) In order to be able to use the LEs H32|l derived 
above, one still needs to calculate the mean collision 
rate C{P)/L, which is determined by Eq. (j7bp : cf. Ap- 
pendix. We also emphasize once again that the approx- 
imation (|25|l . leading to the appearance of C(P), can in 
principle be omitted (in the non-relativistic as well as 
in the relativistic case). More precise results for friction 
coefficients and noise correlations can then be extracted 
from Eq. I|A6|I in the Appendix. 

(iii) The stochastic force x{t) in Eq. (|32|) is (5-correlated 
(memory- free), but non-Gaussian; i.e., in order to com- 
pletely specify the stochastic process one actually has to 
determine all higher order correlation functions. This is 
practically unfeasible. Therefore, in numerical studies 
and/or practical applications, one could use a Gaussian 
approximation (GA) of Eqs. H32I) . obtained in the follow- 
ing manner: We rewrite Eq. Ij32b|l equivalently as 



(33a) 



where 



HP) ^ HP) 



KP) 
p ' 



m^^=^- (33b) 

y^2D{P) ^ ' 



Reminiscent of standardized Gaussian white noise, the 
effective noise force C,{t) is characterized by 



(c(t)) = o, {m~as))^5{t~s), 



(33c) 



but the higher moments are non-Gaussian. Accordingly, 
the GA is achieved by replacing C,{t) in Eq. (|33a|) with 
standardized momentum-independent Gaussian white 
noise C^{t). The resulting stochastic differential equa- 
tion is a standard LE with multiplicative Gaussian white 
noise. Hence, after having specified a discretization rule, 
one can easily write down the corresponding Fokker- 
Planck equation as well as the corresponding stationary 
distribution |6^ . 

The GA obtained this way neglects higher-order cu- 
mulants of the noise, so it cannot be expected that the 
'truncated' LE yields exactly the same relaxation behav- 
ior and/or the same stationary solution as the full rela- 
tivistic LE H32|l ,62b Nevertheless, this approxima- 
tion should provide useful estimates. In particular, if the 
stationary momentum distribution of the Brownian par- 
ticle can be guessed by other arguments |23| , then the GA 
can be made self-consistent with respect to this distribu- 
tion by fixing a suitably generalized Einstein relation for 



the friction and noise coefficients. In this case it suffices 
to calculate, e.g., v{P), because the corresponding noise 
amplitude D{P) is then uniquely determined by the Ein- 
stein relation (i.e., by the stationary distribution). 

Due to the multiplicative noise coupling, the results 
obtained from the GA will also depend on the choice of 
the discretization rule P, 15^, fss!, fsi!, "eS, El El il il . 

Loosely speaking, this discretization dilemma is the price 
that one has to pay for mapping the large number of 
collisions between t and t+r onto a single instant of time. 
Our experience with the non-relativistic LE (cf. remarks 
at the end of Sec. Ill ti\ suggests that the "transport" or 
"kinetic" interpretation, corresponding to the post-point 
discretization rule 54, 62, 63, 64], should be preferable 
in the relativistic case as well. 

(iv) In principle, it should be straightforward to gen- 
eralize the above approach to higher space dimensions, 
by expressing the momentum vector after the collision, 
P, in terms of the momenta before the collision, P and 
p, analogous to Eq. (|21|l . In the 2d or 3d-case, complica- 
tions may arise mostly due to the fact that one also has 
take into account the corresponding collision angles and 
cross-sections (e.g., when determining the collision rates; 
cf. comments at the end of the Appendix). 

(v) According to our above results, the previously pro- 
posed 'relativistic' LEs IMIlllillSfll should be viewed 
as approximations, which can be fruitful for generat- 
ing/simulating ensembles of relativistic particles in a sim- 
ple manner. It is also evident now why these earlier ap- 
proaches were intrinsically limited. Debbasch et al. '3(?| 
have postulated that the relativistic stochastic force in 
the rest-frame of the bath is ordinary Gaussian white 
noise with a constant amplitude Z?, whereas we in our 
prior work 21, 27] assumed Gaussian noise in the (co- 
moving) rest-frame of the Brownian particle. As fol- 
lows from the derivation presented here, neither of these 
assumptions is accurate if one properly takes into ac- 
count both the relativistic conservation laws and the rel- 
ativistic momentum distributions of the heat bath parti- 
cles. However, if we consider suitably chosen momentum- 
dependent friction and diffusion coefficients, then the pre- 
viously proposed 'relativistic' LEs .2L.26, .27. 30] become 
equivalent to the Gaussian approximation of the relativis- 
tic LE 
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APPENDIX A: CALCULATION OF THE 
COLLISION RATE 

We aim to derive an explicit expression for the expecta- 
tion value {Irit,T))^ in the limit r — > 0, as, e.g., required 
in Eqs. (P|). 

By definition, the function Ir{t,T) G {0,1} indicates 
whether or not the Brownian particle has collided with 
the heat bath particle r during the time interval [t, t + t]. 
The positions of the Brownian and heat bath particle at 
time t are denoted by X and Xr, respectively. Ignoring 
the possibility of a collision, for small enough r, the new 
positions at time t + t would be given by 



X' = X + Vt, 



(Al) 



where V and Vr are the velocities. Then, the indicator 
function Ir{t,T) can be explicitly represented as 

Ir{t,T) = e{X ~ Xr) Oix'^ ~ X') e{Vr -V) + 

Q{Xr - X) eiX' - 4) 9(1/ - Vr), (A2) 
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where Q{x) is the Heaviside- function, defined by 

'O, 2; < 0; 
e(x) = <( 1/2, a; = 0; 

1, a; > 0. 



(A3) 



The first (second) summand in Eq. ljA2p refers to the ini- 
tial configuration, where the heat bath particle is located 
at the left (right) side of the Brownian particle. Let us 
list some properties of the collision indicator Ir{t, r). 
First we note that Ir{t,T) is idcnipotcnt, i.e., 

mt,r) ^ Ir{t,T) (A4a) 

holds for 7 = 1,2.. .. Furthermore, for t — > 0, we have 

lr{t,0) = 0. (A4b) 

Accordingly, the Taylor-expansion at r = gives 

'dir 



Ir{t,T) 



dT 



(i,0) 



and some function G{xr,Vr) such that the expectation 
value (G{xr,Vr)'^ exists. We are interested in expecta- 
tions of the form 

as required for calculating the mean value of the stochas- 
tic force and its higher correlation functions [e.g., com- 
pare first line of Eq. (|29() ]. For small r, we may truncate 
the Taylor expansion after the linear term, yielding 

G{Xr,Vr)Ir{t,T))^^(GiXr,Vr)^{t,0)^ T. (A6a) 

Making use of the result l|A4d|l . the mean value on the 
rhs. is given by 



dIr \ 1 f°° 

G{Xr,Vr) -^{t,0)j ^ 2L J '^'"r(Vr-V) X 



(A4c) 



G{X,Vr) fUvr) + 



In order to determine (^(t,0))j^, we note that 



—eix'.^x'] 



r=0 



— e(x, -X + iVr- V)t) 



r=0 



and, analogously. 



^e(x'-< 



{Vr - V) S{Xr -X+ {Vr - V)t) 



^e(X -Xr + iV- Vr)T) 
OT 



{V ~ Vr) S{X -Xr + iV - Vr)T) 
[V -Vr) 6{X ~Xr). 



— J dVr [V - Vr) X 
G{X,Vr) fliVr). 

(A6b) 

In particular, by choosing G{xr,Vr) = 1, wc find the 
=0 collision rate 



r->0 T 



where 



C{V) 



T = 



^it,0)\ =^CiV), (A7a) 



dVr (Vr - V) fl{Vr) + 
V 

AVr {V - Vr) jl{Vr). (A7b) 



Hence, we find 
dl 

-—^{t,0) = iVr-V)e{X -Xr)S{Xr-X)eiVr-V) + 
OT 

(V - Vr) e{Xr - X) d{X ~ Xr) - Vr), 

= 9(0) {Vr - V) S{Xr - X) e{Vr ~ V) + 
8(0) (y - Vr) S{X - Xr) - Vr), 

and, with 0(0) = 1/2, the useful result 
dl 1 

-^(^,0) = -{Vr-V) S{Xr - X) e{Vr ~ V) + 
OT 2 

^{V -Vr) S{X - Xr) Q{V ~Vr). (A4d) 

Now let us consider a spatially homogeneous one- 
particle bath distribution of the form 



The following comments are in order: 

(i) The above derivation is valid for both non- 
relativistic and relativistic heat bath distributions f^{vr)- 
Upon identifying C{P) = C{V{P)), where P is the mo- 
mentum of the Brownian particle, we obtain the rigorous 
justification for Eq. (|7a|l . However, in the non-relativistic 
case we have V — P/M , whereas in the relativistic case 

V = P/{M^ + F2)1/2_ Additionally, we note that the 
support interval of the relativistic velocity distribution 
f^{vr) is given by [— c, c], which determines the effective 
upper and lower integral boundaries in Eq. ljA7b|) . 

(ii) Given a certain bath distribution f^{vr), the ex- 
act result (jA6|l allows for evaluating the quality of the 
product approximations ((HJ and H25(l . respectively. 

(iii) The Stokes approximation corresponds to setting 

V = in Eq. (|A7bp . yielding 



C(0) = - / dVr \Vr\ fliVr). 



(A8) 
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This shows that the Stokes approximation is useful for 
slow Brownian particles, but inappropriate at high ve- 
locities. 

(iv) It is in principle possible to apply the same pro- 
cedure to higher space dimensions, but then the expres- 



sion (jA2|) for the indicator unction has to be modified 
accordingly (e.g., by taking into account the geometric 
shape of the Brownian particle). As a consequence, ana- 
lytic calculations will become much more difficult. 



